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Abstract 

New exact solutions are obtained for several nonlinear physical equa- 
tions, namely the Navier-Stokes and Euler systems, an isentropic com- 
pressible fluid system and a vector nonlinear Schrodinger equation. The 
solution methods make use of the symmetry group of the system in situa- 
tions when the standard Lie method of symmetry reduction is not appli- 
cable. 



1 Introduction 

Lie group theory provides very general and efficient methods for obtaining exact 
analytic solutions of systems of partial differential equations, specially nonlinear 
ones [jlj Q . The different methods have in common that they provide a reduc- 
tion of the original system. This reduction usually means the reduction of the 
number of independent variables occurring, possibly a reduction of the number 
of dependent ones too. 

The "standard", or "classical" reduction method goes back to the original 
work of Sophus Lie and is explained in many modern texts H — H ■ Essentially, it 
amounts to requiring that a solution of the equation should be invariant under 
some subgroup Go C G, where G is the symmetry group of the considered 
system of equations. The subgroup Go must satisfy certain criteria, in order to 
provide such group invariant solutions (see below). 

The purpose of this article is to further develop, compare and apply alter- 
native reduction methods. They have in common the fact that they provide 
solutions not obtainable by Lie's classical method. We shall survey the "tool 
kit" available for obtaining particular solutions of systems of partial differential 
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equations, and further refine some of the tools. In the process we shall obtain 
new solutions of some physically important equations such as the Navier-Stokes 
equations, the Euler equations, the equations describing an isentropic compress- 
ible fluid model and the vector nonlinear Schrodinger equation. 

We shall consider a system of m partial differential equations of order n, in- 
volving p independent variables (x\, x%, x p ) and q dependent variables [u\, «2, u q ) 

A v (x,u ( - n A =0, u = l,...,m, (1) 

where it'™) denotes all partial derivatives of u a , up to order n. 

S. Lie's classical method of symmetry reduction consists of several steps. 

1. Find the local Lie group G of local point transformations taking solutions 
into solutions. Realize its Lie algebra L in terms of vector fields and identify it 
as an abstract Lie algebra. The vector fields will have the form 



p d q d 

v a^^C( x , u )-Q~[ + J2 ( Pa(x,u)—, a = l,..,r = dimG. (2) 

i—l a — 1 

2. Classify the subalgebras L* C L into conjugacy classes under the action of 
the largest group G leaving the system of equations invariant (we have G C G) . 
Choose a representative of each class. 

3. For each representative subalgebra Li find the invariants I 1 {x\, x pi u\, u q ) 
of the action of the group Gi — (exp Li) in the space M ~ X xU of independent 
and dependent variables. Let us assume that p + q— s functionally independent 
invariants I 1 exist, at least locally (s is the dimension of the generic orbits of 
Gi). 

4. Divide (if possible) the invariants I 7 into two sets: {F\, F g }and ...,£ fe }, 
k = p~ s in such a manner that the Jacobian relating and u a (a, fi = 1, q) 
has maximal rank 



9(«i, ...,u q ) ' 



rank J = q. (3) 



Then consider F p to be functions of the other invariants {^jwhich serve as 
the new independent variables and use condition (g) to express the dependent 
variables u a in terms of the invariants. 

5. Substitute the obtained expressions for u a into the original system and 
obtain the reduced system, involving only invariants. The reduced system will 
involve only k — p — s, s > 1 independent variables. 

6. Solve the reduced system. If the variables £* depend on x\,...,x p only, 
this will yield explicit solutions of system (|l|). Otherwise, if depend also on 
the original dependent variables u a , we obtain implicit solutions. 

7. Apply a general symmetry group transformation to these solutions. 
This procedure provides a family of particular solutions that can be used to 

satisfy particular boundary or initial conditions. The classification of subgroups 
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Gi can be viewed as a classification of conditions that can be imposed on the 
obtained solutions. 

Steps 3, 4 and 5 of Lie's method can be reformulated as follows. Take the 
vector fields {v a } forming a basis of the considered Lie subalgebra Qi <Z Q and 
set their characteristics equal to zero 

p 

Qa =l Pa( x : u )-^2Ca(x,u) u^=0, a=l,...,q, a = 1, n = dim 

i=l 

(4) 

Solve the systems (|l|) and (Q) simultaneously. 

Several alternative reduction procedures have been proposed, going beyond 
Lie's classical method and providing further solutions. They all have in common 
that they add some system of equations to the original system (Q) and that the 
extended system is solved simultaneously. These additional equations replace 
the characteristic system (Q). 

In its generality this was proposed as the method of " differential constraints" 
and independently as the method of "side conditions" [jlO|, O. Different 
methods differ by the choice of this system of side conditions. 

Basically two different approaches exist in the literature. The first makes 
further use of the symmetry group G of system (Q), the second approach goes 
beyond this group of point transformations, or even gives up group theory alto- 
gether. Let us briefly discuss some of these methods. 

1. "Group invariant solutions without transversality" . This method was 
proposed by Anderson et al. [Q quite recently and deals with the situation 
when the rank condition (||) is not satisfied. It was shown jj3 that under 
certain conditions on the subgroup Gi C G one can still obtain Gi invariant so- 
lutions. We recall that a group invariant solution, with or without transversality, 
is transformed into itself by the subgroup Gi C G. 

2. The method of partially invariant solutions. A solution u (x) is "partially 
invariant" p, 13 under a subgroup Gi C G of the invariance group, if Gi, when 



acting on u (x) sweeps out a manifold of a dimension that is larger than that of 
the graph of the solution, but less than the dimension of the entire space M. 
A group Gi may provide both invariant and partially invariant solutions (see 
below). However, if the rank condition (^) for the Jacobian is not satisfied, at 
least on a set of solutions, then Gi will not provide invariant solutions, but may 
provide partially invariant ones. If the rank of the Jacobian J is q' , with q' < q, 
then we can express u\, u q i in terms of invariants and let it 9 <+i, u q depend 
on all of the original variables x\,...,x p . We then substitute the dependent 
variables back into the original system and obtain a " partial reduction" . Solving 
this system, we obtain the partially invariant solutions. Irreducible partially 
invariant solutions are partially invariant solutions that cannot be obtained by 
Lie's method using the subgroup Gi or any other subgroup G\ of the symmetry 
group of the considered equations. Such solutions were constructed for certain 
nonlinear Klein-Gordon and Schrodinger equations, in [fl4[ and fl5|| , and for 
some equations of hydrodynamics in Ref . 0, [l?], |l8) . 
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The theory of partially invariant solutions was further developed by Ondich 
[ fl9[ pt| who formulated irreducibility criteria for certain classes of equations. 
For other applications, see Ref. |H| p2| , p3| . 

In terms of the equations (W), the method described above boils down to 
taking only q' < q of the equations (Q). As we shall show, it is possible to 
obtain further partially invariant solutions by different methods. 

Among methods that go beyond the use of the symmetry group G we mention 
the following. 



The Clarkson-Kruskal direct reduction method 24, |25| does not make 



explicit use of group theory. It is postulated that the dependent variables should 
be expressed in terms of new dependent variables that depend on fewer indepen- 
dent ones. The corresponding Ansatz is substituted into the original equation, 
which must then be solved. It has been shown that the direct method is in- 
timately related to the method of conditional symmetrie s p6| , p7| , p8| , to the 
" nonclassical method" proposed by Bluman and Cole |29[ |30|, and to potential 
symmetries [ ^T| . This method can also be interpreted in terms of side condi- 
tions [[l0| . The differential constraints added to system (Q) in this case have 
the form of first order quasilinear partial differential equations of the form (||). 
However the coefficients and £ l a are not related to a Lie point symmetry of 
eq. (|). 

4. The group foliation method. The method goes back to S. Lie, is de- 
scribed by Ovsiannikov [0] and has recently been applied to obtain solutions 
of self dual Einstein equations p2|, |33| . In terms of differential constraints this 
method amounts to embedding the system (Q) into a larger system, consisting 
of all equations up to some definite order, invariant under the same Lie point 
symmetry group as ([!]) (and involving the same variables). 

5. The method of "partial Lie-point symmetries" , proposed by Cicogna and 
Gaeta p4|]. This is a modification of the method of conditional symmetries. The 
method is in some cases easier to use and may provide solutions in cases when 
the equations of the conditional symmetry method prove to be untractable. 

6. For integrable equations |3f| generalized symmetries can be used to gen- 
erate side conditions. These will be higher order equations, rather than first 
order ones. 

7. The method of nonlocal symmetries. This consists in extending the space 
of dependent variables by adding some auxiliary variables, which can be poten- 
tials or pseudopotentials associated to the system of equations under analysis 
[HI H' EOl Hi- -^ nonlocal symmetry will be then a symmetry of the original 
system augmented with the equations defining the new nonlocal variables. 

This article is organized as follows. In Section 2 we introduce the concept 
of "weak transversality" allowing us to simplify the method of "group invariant 
solutions without transversality" |l2| and to relate the rank of the matrix of 
invariants (|J) to that of the coefficients of vector fields. The method of weak 
transversality is applied in Section 3 to obtain new invariant solutions of the 
Navier-Stokes equations and of the isentropic compressible fluid model. In 
Section 4 we establish a relation between partially invariant solutions and the 
transversality condition. This is then applied to obtain new irreducible partially 



4 



invariant solutions of the vector nonlinear Schrodinger equation, of the Euler 
equations, of the Navier-Stokes equations and of the isentropic compressible 
fluid model in (3+1) dimensions. The concept of the irreducibility of partially 
invariant solutions is discussed. Some conclusions are presented in the final 
Section 5. 



2 Group invariant solutions, strong and weak 
trans versality 

S. Lie's classical method of symmetry reduction was outlined in the Introduction. 
The first two steps are entirely algorithmic and we shall assume that they have 
already been performed. Thus, we are given a system of equations ([!]) and have 
found its Lie point symmetry group G, the Lie algebra of which is the symmetry 
algebra L. The symmetry algebra has dimension r and has a basis realized by 
vector fields of the form (0). Each vector field has the property that its n-th 
prolongation annihilates the system ([!]) on its solution set 

prw A u |a m =o= 0, v, (i= l,...,m. (5) 

The functions (x, u) and <^> % a (x, u) are thus explicitly known. 

Let us now consider a subgroup Go C Gand its Lie algebra go- A solution 
u = f (x) of the system ([[J) is Go invariant if its graph Tf ~ {x, f (x)} is a Go 
invariant set: 

S- 17=17, geG . (6) 
The vector field (||) can be written in evolutionary form jlj as 

VE, a = ^2 \P2( x > u )-'52&( x > u ) v %i) a = l,...,r, (7) 

a=l \ i=l / 

A Go-invariant solution will satisfy the qx ro characteristic equations (|J) 
associated with the basis elements of the Lie algebra L$ (dimio = f*o)< 

The following matrices play an essential role in symmetry reduction using 
the symmetry group of the considered system of equations. 

1. The matrices Si and S2 of the coefficients of the vector fields v a spanning 
the algebra L, or its subalgebra Lo, and defined as follows: 

Hi = {£(*,«)}, 5iet rx " 

E 2 = {C(x,u),^(x,u)}, SjeR^' (8) 

where a = 1, ...,r labels the rows, zand a labels the columns. 

2. The matrix of characteristics of the vector fields (||) (or (fjj)) spanning the 
considered algebra L (or its subalgebra Lo) 

Qa = { v aU a } a=l,...,r; a = l,...,q. (9) 
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3. The Jacobian matrix J of the transformation relating the dependent 
variables u a and the invariants of the action of Go on the space M <~ X x U of 
independent and dependent variables. 

Let us now consider a specific subalgebra Lq C L and use it to obtain group 
invariant solutions via symmetry reduction. If the group Go acts regularly and 
transversally on M ~ X x U then 

rank {C a (x,u}} = rank {C a (x, u) , ip% (x, u)} . (10) 

This rank is equal to the dimension of the generic orbits of Go on M. If 
the transversality condition ( |ic| ) is satisfied, at least locally, for all {x,u} € M, 
then Lie's classical reduction method is directly applicable. Indeed, if ( |l0| ) is 
satisfied then the rank of the matrix J of eq. (|3|) is maximal, rank J — q (for a 
proof, see e.g. Ref. jD, Chapter 3.5). It follows that all dependent variables can 
be expressed in terms of invariants and a reduction is immediate (to a system 
with q dependent variables and p — s independent ones). 

If the action of Go on Mis fiber preserving (i.e. the new independent vari- 
ables only depend on the old independent ones), Lie's method provides explicit 
solutions. This happens because the new invariant independent variables z t can 
be chosen to depend only on the original independent variables 

Zi = Zi (xi, ...,x p ) i = l,...,p—s. (11) 

More generally, if we have Zi — Zi (x , u) , we obtain implicit solutions. 

We shall call the rank condition ( |l0| ) " strong transversality" . Quite recently 
[ fl2| a method was proposed for obtaining group invariant solutions when equa- 
tion ( [Toj ) is not satisfied. The method of Ref. |Q can actually be simplified by 
introducing the concept of "weak transversality". 

Definition 1 The local transversality condition will be said to be satisfied in 
the weak sense if it holds only on a subset M C M, rather than on the entire 
space M : 

rank{C a (x,u)} rank {C a {x,u) ,<p% (x,u)} | F/ (12) 

In other words, even if the transversality condition is not in general satisfied, 
there may exist a class S of functions u = f (x) such that for them the condition 
(||) holds. 

The "weak transversality" method is quite simple, when applicable. It con- 
sists of several steps. 

1. Determine the conditions on the functions u — f (x) under which eq. ( [l2| ) 
is satisfied. Solve these conditions to obtain the general form of these functions. 

2. Substitute the obtained expressions into the matrix of characteristics (||) 
and require that the condition rank Q — be satisfied. This further constrains 
the functions / (x). 

3. Substitute the obtained expressions into the system ([!]). By construction, 
the solutions, if they exist, will be Go-invariant. 

This method can only be applied if the matrix elements in the matrix S2 
depend explicitly on the variables u a . This poses strong restrictions on the 
considered algebra go- We shall give some examples of the method in Section 3. 
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3 Examples of invariant solutions obtained by 
the weak transversality method. 

3.1 The Navier-Stokes equations 

The Navier-Stokes equations in (3+1) dimensions describing the flow of an 
incompressible viscuous fluid are: 

Ut + ~u • V Tt + Vp - vV 2 ^u = 0, (13) 

V-~u=0, (14) 

where ~u — (ui (x, y, z, t) , u 2 (x, y, z, t) , u 3 (x, y, z, t)) is the velocity field, p = 
p (x, y, z, t) the fluid pressure and v the viscosity coefficient. 

The symmetry properties of these equations have been intensively investi- 
gated by many authors from different points of viewfsee, for instance, |59| and 



references therein). It is well known |40| that eqs. (13)-(|14|) are invariant under 
the flow generated by the following vector fields: 

(15) 

(16) 
(17) 
(18) 
(19) 

j y -y iu z ~r ^ti^f — u-i^mi ~ "2^u 2 — ^3^3 — 2pd p , (20) 

Li = zd y - yd z + u 3 d U2 - u 2 d u ^ (21) 

L 2 = xd z - zd x + Uid U3 - u 3 d Ul , (22) 

L3 = yd x - xd y + u 2 d Ul - uid U2 , (23) 

where a, /3,7 and d are arbitrary functions of time. The operators ^ generate 
symmetry transformations that can be interpreted as boosts to frames moving 



Bi = 


ad x + 


a d Ul - 




B 2 = 


I3d y + 


P d U2 - 


Pydp, 


B 3 = 


ld z + 


7 d U3 - 


7 zdp, 




T 


= d u 






Q 


= d p , 




zd z - 


- 2td t • 


- uxd Ul 


- u 2 d U2 



with arbitrary velocities v (t) = A ya, P,"f j , where A is a constant. Space trans- 
lations and Galilei boosts are obtained if a, (3 and 7 are linear in t .The operators 
T and Q express the invariance of the eqs. ([l3|)-(|l4]) under translations of time 
and pressure, D generates scaling transformations, and L\.L 2 and L 3 are the 
generators of the group of the rotations of the Euclidean space. 

Example 1. Let us consider the subalgebra generated by Li,L 2 and L 3 . 
Here we apply the ideas discussed in Section 2 to determine rotationally invari- 
ant solutions for the Navier-Stokes equations. 
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The matrices of the coefficients Si — (x, m)) and S2 = Kq u) , 4% {x, u)) 
are represented by 

z -y \ / z -y u 3 -u 2 

—z x I , S2 = — z x — U3 u\ 



y —x I \ y —x —u\ 

(24) 

We observe that the matrix Si has rank 2, whereas the matrix S2 has rank 
3. In this situation, the transversality condition is violated in the strong sense. 
In other words, it is not true that for every function ~Ti = (ui, u 2 , Us) the ranks 
of the matrices Si and S2 coincide. In general, the system of characteristics 
is not compatible, the Jacobian matrix J will not have maximal rank, and it 
is not possible to use the classical symmetry reduction approach. To overcome 
these difficulties, let us force the matrix S2 to be of rank 2. This requirement is 
equivalent to a system of algebraic equations for u , obtained imposing that the 
determinants of all 3 x 3 matrices constructed using the rows and the columns 
of S 2 be equal to zero. Once this algebraic system is solved, we get the class S 
of functions u — (ui, 1^2, U3) on which transversality is weakly restored. The 
class S in this case is characterized by the conditions: 

ui = f (x,y,z,t) x, u 2 = f (x,y,z,t) y, u 3 = f (x,y, z,t) z, p = p(x,y,z,t) 

(25) 

The second step consists in solving the characteristic system Q% (x,uW) = 
for the class S of eqs. (^5|). This forces the functions / and p to have the form: 

f = f(r,t), p = p(r,t), (26) 



where r = \J x 2 + y 2 + z 2 . Relations ( |26| ) represent the most general form for 
the function u and p to be rotationally invariant in the Euclidean space. Sub- 
stituting these expressions for Tt and p into the eqs. (p^)-(|l^), we find the 
solution 

. a (t) a a 2 , 

The same vector fields Li,L 2 and L3 provide also a subalgebra of the sym- 
metry algebra of the Euler equations. Anderson et al.|l2| obtained a class of 
rotationally invariant solutions of the Euler equations by means of their tech- 
nique of reduction diagrams. In Ref. [fl6[ partially invariant solutions related 
to the subalgebra {Li, L 2 , L3} for the equations describing a nonstationary and 
isentropic flow for an ideal and compressible fluid in (3 + 1) dimensions have 
been constructed using the transformation (p5|). A similar situation is also 
observed in magnetohydrodynamics [[l7|. Our solutions of the Navier-Stokes 
equations coincide with the solutions of the Euler equations found in Ref. ]l2| . 
Physically that means that the solutions (^?j) describe a laminar flow for which 
viscosity plays no role. This phenomenon occours because the components of 
the vector Tt in eq. ( p7| ) are all harmonic functions, i.e. they satisfy the Laplace 
equation, in addition to the Navier-Stokes equations. 
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Example 2. Now, Let us analyze the subalgebra denned by the operators 



D = xd x + yd y + zd z + 2td t - u x d ui - u 2 d U2 - u 3 d U3 - 2pd p , 



L3 = yd x - xd y + u 2 d Ul - uid U2 , 
X = t h : d x + kt k - 1 d Ul -k(k-l) t k - 2 xd i 



Y = t k d y + k t k ~ 1 d U2 — k (k- 1) t"-' z yd 1 



fc-2„ 



(28) 



which is a subalgebra of the Galilei-similitude algebra for a given k E R. It is 
immediate to check that the local transversality is violated for this subalgebra 
in the strong sense. The matrix 2 2 is now represented by 



V 



x y 

y ~ x 
t k 
t k 



z 2t 









-Ul 

kt k ~ l 




-U 2 





— U3 — 2p 


-k (k - 1) t--"x 



fc-2 

k(k-l) t k ~ 2 y J 



If we impose that the matrix S2 should have rank 3 (that is weak transversality) , 
we get 

Ui = k — , Ui = k — 
t t 

u 3 = u 3 (x,y,z,t), p = p(x,y,z,t). 

As a second step, let us solve the characteristic system Q" (x, u^ 1 ^) = 0, 
which consists of 16 linear differential equations of first order in the derivatives 
of the velocity components Uj and the pressure p. The most general function 
living in the space of the dependent variables and invariant under the flow 
associated to the generators (p8|) is 



«i 



P = 



u 2 



y 



u 3 



t' f 

k(k-l){x 2 +y 2 ) 

2t 2 z 2 



(29) 



(30) 



with a and (3 arbitrary functions of t z 2 . Substituting into the Navier-Stokes 
equations (|l3])-([l^), we obtain the following three parameter set of solutions: 



Ui = fc-, 1t 2 = k~, 



Ci kz 
u 3 = — = — 2 — 
Vi t 



(31) 

y 2 + Az 2 ) + 2c 2 t} , 
(32) 

with ci, C2, k G M. This solution is invariant under the group generated by the 
algebra (p8|) and satisfies weak (but not strong) transversality. 



P : 



^2 jciVte + k (x 2 + y 2 + AciVtz - 2z 2 ^j - k 2 (a 
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3.2 The isentropic compressible fluid model 



The equations describing the non-stationary isentropic flow of a compressible 
ideal fluid are jy] 

Hi + u ■ VI? + k a Va = (33) 
at + ^u- Va + k~ 1 a V • " 



0. 



(34) 

where u = u\ [x, y, z, t) , u 2 (x, y, z, t) , 1*3 (x, y, z, t) is the velocity field, a = 
a (x, y, z, t) is the velocity of sound, related to the pressure p and the density p 

/ \l/2 

by the formula a = ( ^ J , 7 is the adiabatic exponent and k — 2/(7 — 1). 

The symmetry group G of eqs. (^3|)-(^) was derived in ref. |l8[ |. For k =/= 3, G 
it is generated by the following vector fields: 



Po 

K x = td x 



dt 



d x , P 2 = d y , 



d z 



d Ul , K 2 = td y 



d 



F = xd x 



K 3 = td z 

Li = zd y - yd z + u 3 d U2 - u 2 d U3t 

L 2 = xd z - zd x + uid U3 - u 3 d Ul , 

L3 = yd x - xd y + u 2 d Ul - uxd U2 , 

yd y + zd z + td t , G = -td t + uid Ul + u 2 d U2 4 

We mention that for k = 3 the symmetry algebra contains an additional 
element generating projective transformations. The operators Pi , Ki and Li are 
the infinitesimal generators of space translations, Galilei boosts and rotations, 
respectively. The operators F and G generate scaling transformations. 

Example 3. Let us consider the subalgebra {L 3 , F + G, K\, K 2 } ■ The ma- 
trix S2 is given by 



u 3 d., 



(35) 

(36) 
(37) 
(38) 
(39) 

„, (40) 



y -x 
x y 



t 

V 



u 2 

z u\ 

1 





-Ui 



"2 



1 



\ 

u 3 a 



/ 



(41) 



and the transversality is again violated in the strong sense, because rank Si = 3 
and rank S2 = 4. If we force the matrix S2 to be of rank 3, then we get the 
following constraints 



x 



y 

U 2 = J, 



u 3 



u 3 (x,y,z,t), a = a(x,y,z,t) . 



From the characteristic system (2, uW) = we deduce 
u 3 = zW(t), a= zA (t) , 



(42) 



(43) 
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where W and A are arbitrary functions of time. Now, substituting relations 



(|2h-(43) into the system (B3h-(M) we obtain the relation 



= \!-\i W2 + W ') ( 44 ) 



and a second order ODE for W 
W 



" + 2 ^2 + WW' + 2 ^1 + -H W 3 + (W + W 2 ) = (45) 

In general, eq. (EJ|) does not have the Painleve property. For special values 
of the parameter k, namely k = —1 and k = —2, it does. In these cases it can 
be reduced to a canonical form (see Ref. p2[ , p. 334) via a linear transformation 
of the type 

W = a{t) U(z(t))+0(t). (46) 
1. For k = — 1 we have 

W" = —2WW' + p (t) (W' + W 2 ) , (47) 

where p (t) = A/t. This equation can be integrated and its solution, which is 
regular, is 

W = V\ 7 \ — » (48) 

where ci and C2 are constants and I n (x) is the modified Bessel function of the 
first kind. Correspondingly we find 



A = Cl t 2 (49) 



in eq. (Q). 

2. For k = —2 we have 



W" = -ZWW - W 3 + q (t) (W' + W 2 ) , (50) 
where q it) = 2/t. In this case, we can integrate and the general solution is: 



W= /* 3 + Cl . (51) 



The solution of eq. (Hi) is 



A = 2V3\-. , (52) 

where ci, and C2 are constants. The solutions for k = — 1 and k = —2 represent 
nonscattering waves. 
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4 Partially invariant solutions of systems of dif- 
ferential equations and the transversality con- 
dition 

A useful tool applicable to the study of systems of differential equations, and 
intimately related to the standard Lie approach, is the theory of partially in- 
variant solutions. The relevant notion in this context is the defect <5of a k— 
dimensional manifold M with respect to a Lie group G. When the group acts 
on a p-dimensional submanifold T C M, it sweeps out an orbit G(T). The 
manifold T will be identified with the graph of a function u — f (x), so its 
dimension will coincide with the number of independent variables (also denoted 
p). As we already said, the case G (Tf) — Tf corresponds to the G-invariance 
of the manifold. Otherwise, G(r/)will be a more generic subset of M. There 
is no guarantee that this subset will be a submanifold. However, if the inter- 
section between an orbit O of G and Tf has a dimension which is constant in 
a neighbourhood N of a point of Tf, then there exists a neighbourhood G of 
the identity of G such that the subset G (N n T f) is a submanifold Q. In the 
subsequent considerations, G (Tf) will be considered as a submanifold. 

Let G be a group, acting regularly with s-dimensional orbits. We call the 
number 

5 = dhaG(T f )-dimTf (53) 



the defect 8 of the function / with respect to G . The usual G-invariant functions 
correspond to the case S — 0. A function will be said to be generic if 5 = m = 
min {s, k — p} . The more interesting situation is when < 6 < too, which is the 
case we will dealing with. In this case, the function / will be said to be partially 
invariant ||. 

Let us consider the system (^) of partial differential equations, whose sym- 
metry group G acts on the p + g-dimensional space M = X X U. Let g be 
a subalgebra of the symmetry algebra of A, and Q the characteristic matrix 
associated to the set of its generators. Then u = f (x) is a partially invariant 
solution of A with defect S with respect to 9 if and only if || 

rank (q (x,it (1) )) = 5. (54) 

The condition ( [si]) provides a system of differential equations involving the 
dependent variables u = (ui,...u q ). In order to determine partially invariant 
solutions, we can extend the original system A by adding the set of differential 
constraints given by the condition (|54|). We must then solve the extended system 
consisting of eq. (Q) and (Q). The set of equations given by the prescription 
( p4[ ) is less constraining than the set required to obtain G-invariance, as in 
formulas (||). 

In this section we will study the role of the local transversality condition 
(and in particular of the notion of weak transversality) in the theory of partially 
invariant solutions and propose a strategy to find them. 
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Let us start by noticing that a partially invariant solution of a system of 
differential equations can be naturally related to the violation of the transver- 
sality condition. Indeed, let A (2, u^ n A = be a system of differential equations 
defined over M C X x U and u = u (x )be a solution of A. Let G be an 
r-dimensional subgroup of the symmetry group of A, acting regularly on M, 
whose generators are given by (Q). If the condition 

rank (C a {xo, uo)) < rank (f* (so, «o) , C (xo,u )) (55) 

is satisfied, then u = f (xq) is a partially invariant solution of A (or possibly 
a generic one). 

Example 4. The vector nonlinear Schrodinger equation 

iip t + Aijj = (^V) ip, (56) 

where i/j 6 C N and A is the Laplace operator in n dimensions, plays an 
important role in many areas of physics. For instance, in nonlinear optics it 
describes the interaction of electromagnetic waves propagating with different 
polarizations in nonlinear media |45|| . In hydrodynamics, it furnishes a model 
for the description of the interactions of N water waves in a deep fluid |46|— f48|. 
For these and other applications of the vector nonlinear Schrodinger equation, 
see also Ref. g§. 

Let us consider the case of three components (N — 3) and two spatial dimen- 
sions (n = 2). The symmetry algebra has been computed in Ref. J50| | . In terms 
of amplitude and phase, the components of the wave function will be written as 

— p i & u}i ■ In particular, we will discuss the role of the subalgebra 

{d x , dy, yd x - xd y + aidui + a 2 duj 2 + 03^3} , (57) 

generated by the two translations in the plane and a rotation combined with a 
transformation of the phases. The matrix of the coefficients H2, given by 

1 
1 

y —x a\ 02 03 

has rank 3, unless (a%, a%, 03) = (0,0,0). Let us consider the case a 2 = 03 = 0, 
a% 7^ 0. We obtain p i — p i (t), (i = 1,2,3), w\ = u>i (t) ,lu 2 = u>2 {t), but U3 is 
not an invariant, so we keep UI3 = 0J3 (x, y, z, t). Substituting into eq. (|56|), we 
obtain: 

Pl = 7tt= i ? ' ^2 = 7 2 > Ps = 7 3 > (58) 

^2 = ^3 = ^ In— --(7I+73)*- (60) 
to t — to 
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This solution is partially invariant with respect to the subgroup corresponding 
to the subalgebra ( |57| ) with a 2 = (13 = and a\ ^ 0, and is not reducible (unless 
we choose to = 0). For to = it is invariant under rotations in the xy-plane. 
Example 5. Let us consider again the isentropic compressible model (for- 



mulas (|33|)-(p 
given by 



3)). For the subalgebra g = {K±, K2, -K3, P3} the matrix H2 



is 



\ 



Here transversality is violated also in the weak sense, 
corresponding Lie group are F = u\t — x, G = u 2 t — y, 
of (|3|) is not invertible, but we can write 



The invariants of the 
a and t. The matrix J 



ui 



F(t) 



t 



U2 



G(t)+y 



A(t), 



(61) 



but leave u 3 general, i.e. u 3 = 
341) we obtain the solution fllf 



U3 (x, y, z, t) . Substituting eqs. (61) into ( 



Ml 



x 



u 2 



u 3 



t + t Q 



t 2 (t + 1 ) 



(62) 



where ^ = j , ^ 2 = f 1 ^ i s an arbitrary function of ^ and £ 2 j c and to are 
constants. The rank of the matrix Q™ of the characteristics associated to ( |62| ) 
is equal to one, and therefore this solution is partially invariant with respect 
to the subalgebra g, with 5=1. Now let us check if it is reducible under any 
other subalgebra of the full symmetry algebra (|35|)-(f40|). To do this, it is useful 
to study the kernel K of the characteristic matrix Q for the full symmetry al- 
gebra (p5|)-(^0|) associated to the solution j6^ ) and to determine its generators 
{ki, k/}. It is clear that if at least a subspace of K can be generated by con- 
stant vectors, then the solution will be reducible with respect to the subalgebra 
identified by these vectors. 

In the case of the solution (^32|), the kernel is generated by 8 vectors, each 
having 12 components. It is possible to show that there exists only one constant 
generator, namely 

k = (0,0,0, t , 0,0, 0,0, 0,1,0,0). 

This implies that the solution ( |62| ) is reducible with respect to the one-dimensional 
subalgebra {K 3 + toP^}- 

The concept of " irreducibility" of a partially invariant solutions needs further 
clarification. Once a solution u = / (x), partially invariant under Gi is found, it 
is of course possible to verify whether it is invariant under some other subgroup 
G'i C G. Let this be the case, let G\ satisfy the strong transversality condition, 
and have generic orbits of dimension r[. The standard Lie method using the 
subgroup G- would then reduce system (|l|) to a system with q — r[ variables. 
In general, specially for q — r[ > 1, we may not be able to solve this system 
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and the methods of partial invariance for the original subgroup Gi may be more 
tractable. If the invariance subgroup G\ does not satisfy the weak transversality 
condition, it may not help us at all. 

We observe that the violation of the transversality condition is not a nec- 
essary condition for the existence of partially invariant solutions of a system 
A (x, u^ 71 ') = 0. In fact, there could be solutions of A which are not solutions 
of the characteristic system , even if it is compatible as an algebraic system. 

Example 6. As a matter of fact, a counterexample was provided by Ondich 
[E9L namely the two variable Laplace equation, expressed in the following form: 



(63) 



This system is clearly invariant with respect to the translations in the plane, 
generated by the vector fields {d x ,d y }. The characteristic matrix Q has the 
form 

-u x -v x -w x \ (g4) 

~Uy ~Vy ~Wy J 

Invariant solutions are obtained if and only if the rank of this matrix is equal 
to zero. This implies that u, v and w are constants: 

u = k, v = A, w = /i, k, A, /x € M. 

Let us impose that the rank of Q is equal to 1. This means that the two 
rows in eq. ( |64| ) are proportional. Then solving the corresponding equations 
and replacing the result into eqs. ([33]), we get the following solution: 

(u, v, w) = (a x + b y + c, a, 6) , a, b E K. (65) 

By construction, this solution is partially invariant with respect to the group 
of translations in the plane, with defect S = 1. Nevertheless, the transversality 
condition is satisfied. More generally, the transversality condition is always 
satisfied if we have rank (£ % a (x, u)) = dim L. 

We mention that the previous solution (|6^) can be also obtained starting 
from the symmetry subalgebra {d x , d y , d u } for which both the weak and strong 
transversality conditions are violated. In this case, the invariants are Ii = v = a, 
1-2 = w = b, where a and b are arbitrary constants. The function u remains 
arbitrary. Putting these constraints into eqs. (|63|), we recover immediately the 
solution (|65|). However, it should be noticed that the analysis of the rank of the 
characteristic matrix furnishes a complete characterization of the invariant and 
partially invariant solutions of a system of differential equations, and provides a 
more general procedure than the use of the invariants. For instance, this second 
approach does not allow to recognize the solution ( |65[ ) as a partially invariant 
one with respect to the subalgebra {d x ,d y }. 

Example 7. To show that a symmetry group satisfying the strong transver- 
sality can produce both invariant and partially invariant solutions, let us con- 
sider the Euler equations for an incompressible nonviscuous fluid in (3+1) di- 
mensions: 

ut + It ■ VI? + Vp = 0, (66) 
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The symmetry group of the Euler equations is well-known jl| ^. It coincides 
with the symmetry group of the Navier-Stokes equations, except that it contains 
an additional dilation. Thus, D of eq. ( pp| ) is replaced by 

Di = xd x + yd y + zd z + td u (68) 

D 2 =td t - uid Ul - u 2 d U2 ~ u 3 d Ua - 2pd p . (69) 

We consider here the subgroup of Galilei transformations. Its Lie algebra is 
given by 

K x =td x + d Ul , K 2 =td y + d U2 , K 3 =td z + d U3 . (70) 

Here the transversality holds in the strong sense and indeed an invariant solution 
will have the form 

ui = ~+Fi(t), U2 = !+F 2 (t), u 3 = j + F 3 (t), p = P{t). (71) 



Let us now look for partially invariant solutions of the system (|5^-(|67j) with 
respect to the same subgroup and impose that the defect be 5 — 2. Writing 
down the characteristic system associated to ( f70| ) and imposing that rankQ = 2, 
we get the following constraints on U\ and u 2 

x z 

U\ = fi\— + pt\u 3 + hi (t) (72) 

u 2 = fiu 3 + -- fx- + h 2 (t) (73) 

but u 3 = u 3 {x, y, z, t) and p = p (x, y, z, t) remain arbitrary. 

Substituting the relations (|72|)-([73|) into the Euler equations, and choosing 
for simplicity h\ — h 2 = 0,wc obtain the solution: 

Ul = t [vfi (1 + A 2 ) + 1] {X ^ (1 ~ 2A2) + 1] ~ 3AM (W + + 

9 ^ f \y — x \ u z — x\ ,„,, 

A ^(-V--V-J' (74) 

U2 = t hfl (1 + A 2 ) + 1] {V ^ (A ' " 2) + 1] " 3M ( VX + Z)} + 

4.2 TP ( X V ~ X ^^Z-X\ 

—J' (75) 



u 3 = rr^^ {« [M 2 (l + A 2 )-2]-3(A, + , )} + 
9 ^ ( Av — x \uz — x\ 

t F [ JL l—>- !L -L > ( 76 ) 
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^- >(aV + m' + i) ( Ai + ' + £> +v{,) - (77) 

where A, /x are constants, i* 1 is an arbitrary function of its arguments and p is an 
arbitrary function of t. We have checked explicitly that, if F is kept arbitrary, 
this solution is not invariant under any subgroup of the symmetry group. Thus, 
it represents an irreducible partially invariant solution of the Euler equations 
of defect 6 = 2, with respect to a subgroup satisfying the strong transversality 
condition. If 

where $is an arbitrary function of its argument, then the solution (|74|)-(f77|) is 
invariant under the subgroup generated by aD\ + 6D2, where D\ and D2 are 
defined in eqs. H])-®. 

However, this subgroup provides a reduced system with three independent 
variables that would be very difficult to solve. 

Particularly interesting is the case when partially invariant solutions can be 
found that satisfy weak but not strong transversality. Indeed, imposing weak 
transversality basically means that a class of functions u = f (x) is chosen 
in such a way that the characteristic system is algebraically compatible. This 
condition is of course not sufficient to guarantee the invariance of these functions 
under the action of the considered group G. Indeed, if we compute the rank of 
the matrix on this class of functions, in general it will be not equal to zero. 
Then, once weak transversality is satisfied, we can choose either to have group 
invariant solutions, using the method outlined in Section 3, or to use the class 
of functions u = f (x) to get partially invariant solutions. 

In the next two examples, we will see how this approach can be used to obtain 
in a simple and straightforward way new classes of solutions of hydrodynamic 
systems. 

Example 8. In Example 2 we studied the algebra ( p8j ) which is a subalgebra 
of the symmetry algebra of both the Euler and the Navier-Stokes equations. We 
showed that the requirement of weak transversality implies 

ui = fcy, it 2 = k~, (78) 

U3=u 3 (x,y,z,t), p = p(x,y,z,t). (79) 

These formulas define a class of functions which is partially invariant with defect 
6 = 2 with respect to the subalgebra (^8|). At this stage, we can choose either to 
have group invariant solutions or partially invariant ones. Indeed, in Example 
2 we forced the class of functions (|78])-(|79]) to be a solution of the characteristic 
system, and then substituting the obtained expressions (|29|)-(|30|) into the system 



13|)-(14) we constructed the group invariant solutions (|3 1|) — (J32 

Another possibility is to substitute formulas (f78|) -(|79"| ) directly into the Euler 
equations (|66|)-(|67|) or the Navier-Stokes equations (|13|)-(|l4|) without requiring 
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further invariance properties. For the case of the Euler equations we get the 
constraints 

t 2 p x + k (k - 1) x = 0, (80) 
t 2 P y + k{k-l)y = Q, 



u 3z 



2k 

T = ' 



u 3 t + u 3 u 3z + - (x u 3x + y u 3y ) +p z = Q. 



(81) 
(82) 

(83) 



We solve this system and obtain the following partially invariant solution of the 
Euler equations: 

ui = k—, U2 — k—, (84) 



t 



t 



2kz n „ ( i y \ 

u 3 = — + x 2 F 

t V xJ 

k (k- 1) (x 2 +y 2 ) k(2k + l)z 2 



(t), 



(85) 



(86) 



2t 2 t 2 

where £ = tx~i,rj — ^ and / (t) and F (£,77) are arbitrary functions of their 
arguments. This solution is irreducible for a generic function F. If F satisfies 
the equation 



ci - 



k 



+ c 2 £ - c 3 



in _ C4 k+1 



c 5 e 



Ft 



k k 

-cs (if + 1) - C4t k V + c 5 e fe J F„ + (2 Cl + c 2 + 20377) F — c 6 (4fe + 1) £ 2fc = 0, 
then it is invariant under the subgroup generated by 

X = ci-Di + c 2 L>2 + C3L3 + c 4 Bi + c 5 S 2 + c 3 B 3 , 

where ci,...,cq are real constants, D% and D 2 are the dilations given by eqs. 
(38)— (69), L 3 is the generator (|2^), and the functions appearing in the boosts 
(15)— (17) are now monomials in t, namely 



a = f 



P = t , 7 = t 



2k+l 



The same procedure can be applied to the Navier-Stokes equations. Repeat- 
ing the previous steps, we obtain the following solution: 



u\ = k—, 
t 



2kz 



"3 



i,y 

u 2 = k-, 



+ a(x,y,t) 



k (k - 1) (x 2 + y 2 ) k(2k + l)z 2 



2t 2 



t 2 



/(*), 



(87) 
(88) 
(89) 
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where a (x,y,t) satisfies the following equation 



k 

a t + -{xa x +ya y -2a)-v (a xx + a yy ) = 0. 



(90) 



We thus obtain a large class of partially invariant solutions of the Navier-Stokes 
equations parametrized by the solutions of the linear partial differential equation 
©■ 

Example 9. Partially invariant solutions with weak trans versality can 
be also found for the case of the compressible fluid model (33)-(|3~i|). Let us 
again consider the subalgebra {L 3 , F + G, K\, K2} and the corresponding weak 
transversality condition (42). Substituting into equations (^)-(Q) we obtain 
a coupled system of quasilinear first order PDE's, namely: 



0. 



0, 



"3i + u 3 u 3z + -u 3x 



a t +u 3 a z + 



t 



-u 3y 



u 3 . 



kaa~ 



0. 



(91) 



In particular, if we assume a z = 0, we reobtain the solution (p2\). However, 
the system (|9l|) allows much more general solutions. 



5 Conclusions 

The main conclusion of this article is that one can do considerably more with 
the symmetry group G of a system of partial differential equations than apply 
the standard method for finding group invariant solutions. 

Indeed, let us assume that the largest group G (of local Lie point trans- 
formations) leaving the system (|l]) invariant has been found and its subgroups 
classified. For each subgroup Go, or its Lie algebra Lq we should proceed as 
follows. 

1. Obtain invariant solutions. First check whether the transversality con- 
dition ( pL| ) is satisfied (in the strong sense). If it is, we apply Lie's classical 
method. This is always possible since transversality assures that the rank con- 
dition ([|) is satisfied. If the (strong) transversality condition is not satisfied, we 
may still be able to obtain solutions invariant under Go |ll| , by imposing " weak 
transversality" on solutions, as described in Section 2 above. This is illustrated 
in Section 3 by Examples 1,2 and 3. 

2. Obtain partially invariant solutions. These can be obtained by (at least) 
three complementary methods. If the transversality condition ( |l0|) is not satis- 
fied and weak transversality cannot be imposed, then the characteristic system 

is not consistent and the rank condition (0) for the invariants is not satis- 
fied. We then choose a subset of q' < q invariants such that we can express q' 
dependent variables in terms of invariants. The remaining q — q' variables u a 



19 



are considered as functions of all the original variables x%,...,x p . Examples 4 
and 5 of Section 4 are of this type, as are those of Ref . Q , Q . 

If the transversality condition is satisfied for Go we may still be able to 
obtain partially invariant solutions, in addition to the invariant ones. Instead 
of imposing that the matrix of characteristics (Q) have rank zero (i.e. that all 
equations (Q) be satisfied) we require that on solutions we have rankQ = 8, with 
8 = 1,2, as the case may be. This rank condition must be solved explicitly 
for u a and the result substituted into eq. (0) (see Examples 6,7). The third 
possibility is to impose weak transversality (if possible), but still to impose 
rankQ = 8 > 1 (see Examples 8 and 9). 

3. Go beyond invariant and partially invariant solutions, either by the 
method of group foliation ||, ||[ or by other methods not using the sym- 
metry group G H-Q. 

Missing at this stage are clear criteria that tell us which approach will be 
fruitful. Furthermore, the same solution may be obtained by different methods 
and it is not clear which of these method will lead to the least amount of 
calculations. For instance, solutions partially invariant under some subgroup 
d C G have been called "reducible" ||-|o|, if they are actually 

invariant under some other subgroup G 2 C G. However, it may be more difficult 
to use G2 , than G\ , specially if the dimension of G2 is small with respect to the 
number of independent variables. 
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